New quantum phase transitions in the two-dimensional Ji — J2 model 
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We analyze the phase diagram of the frustrated Heisenberg antiferromagnet, the Ji — J2 model, 
in two dimensions. Two quantum phase transitions in the model are already known: the second 
order transition from the Neel state to the spin liquid state at (J2/Ji)c2 = 0.38, and the first 
order transition from the spin liquid state to the coUinear state at [J2/ Ji)c4, = 0.60. We have found 
evidence for two new second order phase transitions: the transition from the spin columnar dimerized 
state to the state with plaquette type modulation at {J2/ Ji)cS. = 0.50±0.02, and the transition from 
the simple Neel state to the Neel state with spin columnar dimerization at {J2/ Ji)c\ = 0.34 ± 0.04. 
We also present an independent calculation of (J2/Ji)c2 = 0.38 using a new approach. 



PACS: 75.10.Jm, 75.30.Kz, 75.40.Gb, 75.30.Ds 

The nature of the quantum disordered phases of low- 
dimensional quantum antiferromagnets is a topic of fun- 
damental importance for the physics of quantum mag- 
netism Such phases can result from mobile holes in 
an antiferromagnetic background as in the t — J ov Hub- 
bard model at finite doping. Alternatively, competition 
of purely magnetic interactions can also lead to destruc- 
tion of long-range order. A typical example of the second 
kind is the Ji — J2 model which exhibits a quantum dis- 
ordered (spin-liquid) phase due to second-neighbor frus- 
trating interactions. Even though it has been intensively 
studied during the last ten years, the J\ — Ji model appar- 
ently still holds many secrets. This model is also an ideal 
testing ground for the theory of quantum phase transi- 
tions because it has very complex dynamics and contains 
a variety of transitions. Exact diagonalization studies ||] 
have shown that the excitation spectrum of the model 
is quite complex and that finite-size effects are large [1). 
Spin-wave like expansions around the simple Neel state 
(which occurs for small frustration) naturally cannot give 
any information about the ground state at stronger frus- 
tration, and consequently non-perturbative methods are 
needed to analyze the latter regime. 

An important insight into the disordered regime was 
achieved by field-theory methods and dimer series 
expansions The above works have established 

the range of the disordered regime, 0.38 < g < 0.60 
(.9 — J2/J1), and have also shown that the ground state 
in this regime is dominated by short-range singlet (dimer) 
formation in a given pattern (see Fig.l). The stability of 
such a configuration implies that the lattice symmetry 
is spontaneously broken and the ground state is four- 
fold degenerate. This picture is somewhat similar to 
the one dimensional situation, where the Lieb-Schultz- 
Mattis theorem guarantees that a gapped phase always 
breaks the translational symmetry and is doubly degen- 
erate, whereas gapless excitations correspond to a unique 
ground state 

Two very recent calculations [§,0 performed by Green 
function Monte Carlo methods have raised new questions 
on the structure of the intermediate phase. The authors 



of Ref. p claim stability of the "plaquette RVB" state 
at g « 0.5. Reference comes to a different con- 
clusion: there is a columnar spin dimerized state with 
plaquette type modulation along the columns. An addi- 
tional very interesting observation ]l0| is that the colum- 
nar spin dimerization penetrates into the Neel phase to 
g « 0.3. To conclude the list of observations which do 
not agree with a simple spin liquid with columnar dimer- 
ization we mention the divergence in the plaquette sus- 
ceptibility found in Ref. Q] at 5 ~ 0.5. 

In the present paper we elucidate all the above ques- 
tions and come to the conclusion that two additional 
quantum critical points exist in the phase diagram of 
the system. These critical points correspond to a new 
generic type of second order quantum phase transition 
considered in Ref. jl^]. At each of the critical points 
points there is condensation of some singlet excitation 
and the critical dynamics is described by the nonlinear 
0(1) cr-model. 

The Hamiltonian of the Ji — J2 model reads: 



H 



Ji Si ■ Sj + J2 Si ■ Sj 



(1) 



where Ji is the nearest-neighbor, and J2 is the frustrating 
next-nearest-neighbor Heisenberg exchange on a square 
lattice (see Fig.l). Both couplings are antiferromagnetic, 
i.e. Ji,2 > and the spins Si = 1/2. We also use the 
notation g — J2/J1. The spin columnar dimerization at 
9 > gc2 is well established ^-0] and therefore we start 
our consideration from this state shown schematically in 
Fig.l. If there is an instability with respect to some kind 
of additional ordering then the gap in the spectrum of 
some singlet excitation must vanish at the corresponding 
critical point jl^. We do not have a reliable technique 
for direct calculation of the singlet gap, but we do have a 
well developed series expansion technique for calculation 
of static susceptibilities. A static susceptibility is pro- 
portional to the corresponding Green's function at zero 
frequency 



Xq oc Gq(cj = 0) 



(2) 



1 



where ujq is the quasiparticle energy, and Zq is the quasi- 
particle residue. So at the critical point 1/x must vanish 
approximately as {g — gcV , with 7 ~ 2{v — rj), where ly 
is the critical index for the spectral gap, A (x (g — gcY ^ 
and 27] is the critical index for the quasiparticle residue, 
Z^{g~ gcT^ 

To analyze possible plaquette type modulation we cal- 
culate the susceptibility of the spin columnar dimerized 
state with respect to the field 



(3) 



2J 



which breaks the translational symmetry in the direction 
perpendicular to the dimers. The series has been com- 
puted up to the seventh order in the dimerization param- 
eter. Results for 1/xp are shown in Fig. 2. The value of 
1/xp vanishes at gcs = 0.50 ± 0.02 and this is the criti- 
cal point for the second order quantum phase transition 
from a simple columnar dimerized state to the 8-fold de- 
generate columnar dimerized state with plaquette type 
bond modulation in the direction perpendicular to the 
dimers suggested in Ref . . This phase transition is of 
the generic type considered in Ref. and therefore it 
is described by 2D nonlinear 0(1) a-model. The critical 
indexes for this model are [Q: v « 0.630, rj « 0.034. 
Therefore one shall expect 7 — 2{v — rj) ^ 1.19. On the 
other hand the Dlog Pade approximants to the series xp 
give 7 = 0.9 ± 0.1. This is fair agreement, and we of- 
fer an explanation for the small discrepancy. The phase 
transition is related to the condensation of some singlet 
excitation which can be considered as a bound state of 
triplet excitations. 



\s) = a2\tt) + a^lttt) + a4\tttt) + .. 



(4) 



We would like to stress that there is very strong mixing 
between two-triplet and multi-triplet bound states. This 
mixing was the reason why vanishing of the singlet gap at 
9 — ffc3 was missed in Reference ^ . In Ref. analysis of 
the singlet excitation was based on a two-particle Bethe- 
Salpeter equation with further account of multi-particle 
contributions as a small perturbation. This assumption 
was wrong because of the strong mixing. So at 5 = gc3 we 
have condensation of effectively a multi-particle bound 
state with relatively small two particle component. The 
mixing between two-particle and multi-particle compo- 
nents of the singlet excitation varies with g and this effect 
cannot be taken into account in the non-linear cr-model 
which assumes condensation of an "elementary" (=struc- 
tureless) field. Ultimately very close to the critical point 
the variation of the mixing can be neglected and one shall 
expect restoration of the pure cr-model behavior. How- 
ever it happens in so narrow vicinity of the critical point 
that the present numerical data cannot assess it. 

Let us consider now the appearance of spin dimer order 
at g = gd as g is increased from small values. A scenario 



put forward some time ago Q and based on the analy- 
sis of the Sp{N), N ^ 00 field theory suggests that the 
dimer order appears simultaneously with disappearance 
of the Neel order, gd = gc2- The dynamics in the vicin- 
ity of the critical point is described by the nonlinear 0(3) 
cr-model in spite of an additional dimer order parameter. 
The additional gapless excitation is irrelevant to the crit- 
ical dynamics because this excitation has extremely large 
size [Q: r ~ l/(g — 5c)*^, M ^ 1. Another possibility is 
that gd < 5c2 and hence there are two separated quan- 
tum phase transitions p3) . The transition at gc2 = 0.38 
is still described by the nonlinear 0(3) cr-model, while 
the transition at gd is of the 0(l)xO(l)-type. So in 
the vicinity of the point gd there is an additional effec- 
tive singlet field which can condense either at momentum 
k = (ti'jO) or (0,7r). The sign degeneracy of the scalar 
field together with the momentum degeneracy gives a 
four-fold degenerate ground state which exactly corre- 
sponds to the degeneracy of the spin-dimerized state. A 
recent work based on the Green function Monte Carlo 
method |l§ gives a hint in favor of this picture. 

Let us give the precise meaning to the terms relevant 
and irrelevant singlet excitation. We consider a quantum 
critical point at which the singlet gap As vanishes. An 
external field which is coupled to the singlet excitation, 
(s|i^|0) 7^ 0, is applied to the system. If the correspond- 
ing susceptibility given by eq. (^) is diverging at the 
critical point we call this singlet excitation "relevant" . If 
the susceptibility is not diverging we call the singlet ex- 
citation "irrelevant" . It is clear that for an irrelevant ex- 
citation the quasiparticle residue Z vanishes faster than 

To analyze the problem of spin dimer order we calcu- 
late the susceptibility of the Neel state with respect to the 
external field which probes spin columnar dimerization. 
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In this calculation we use the usual Ising series expansion 
up to seventh order. Note that in spite of the similarity 
between ^ and (|5) these are two quite different situ- 
ations. The field (|) assumes that the dimers aligned 
along the i-direction already exist and it probes a possi- 
ble modulation in the j-direction. The field (|^) is applied 
to the Neel state and therefore it does not assume any 
dimer order. The values of 1/xd versus g are plotted in 
Fig. 3. It is clear that 1/xd vanishes somewhere in the 
interval g G [0.3, 0.4], but the data is not precise enough 
to distinguish gd from gc2- To distinguish between the 
two scenarios discussed above we have to realize that Fig. 
3 clearly indicates the relevant singlet excitation. In the 
case of an irrelevant singlet Q the quasiparticle residue 
is extremely small, Z oc {g — gc)^ , M ^ 1 and hence 
the susceptibility has no divergence at the critical point. 
Thus we conclude from Fig. 3 that gd = 0.34 ± 0.04, 
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and that gd < gc2, so there is a region Qci < 9 < 5c2 
where the spin columnar dimer order and the Neel order 
coexist. The critical dynamics at Qci is described by the 
relevant gapless singlet excitation. There is no doubt that 
the irrelevant gapless singlet excitation at g w gc2 also 
exists, but it has an exponentially small residue and 
hence its contribution to the susceptibility is negligible. 

The final result we report here is a new way of esti- 
mating gc2 — 0.38. The previous best calculation j|] was 
based on vanishing of the triplet gap in the spin liquid 
phase. A previous attempt to estimate gc2 by Ising 
expansions for the staggered magnetization in the Neel 
phase showed the magnetization vanishing around 0.4, 
but the series were erratic in this region and the preci- 
sion low. The new estimate is based on Ising expansions 
Ip^ in the Neel phase for the 1st, 2nd and 3rd neighbor 
correlators (SfSj) and {SfS^), where z is the direction of 
staggered magnetization. The series has been computed 
up to order 9 for 1st and 2nd neighbor correlators and 
to order 7 for 3rd neighbor correlator. The differences 
of these correlators are shown in Fig. 4. The transition 
point is identified by the condition {S'fS^) = (S^Sj), 
corresponding to restoration of spin rotational symmetry 
to the ground state. This gives gc2 — 0.38(3), in excellent 
agreement with previous results. 

In conclusion, the zero temperature phase diagram 
and the excitation spectra of the Ji — J2 model are 
shown schematically in Fig. 5. There are four critical 
points: gel = 0.34 ± 0.04, 5^2 = 0.38, 3^3 = 0.50 ± 0.02, 
gc4 — 0.60. The states are: g < gd - the simple Neel 
state, gd < g < gc2 - the columnar dimerized Neel state, 
5c2 < 9 < 9c3 - the simple columnar dimerized spin liq- 
uid, gc3 < g < gc4 - the columnar dimerized spin liquid 
with plaquette type modulation, g > gd - the coUinear 
state. The transitions at gd and gcz are second order 
phase transitions of the 0(1) x 0(1) and 0(1) symmetry 
classes correspondingly. Energies of the relevant singlet 
excitations vanish at the critical points. The transition 
at gc2 is a second order phase transition described by the 
nonlinear 0(3) cr-model. The energy of the triplet exci- 
tation vanishes aX g < gc2', at the critical point there is 
also a singlet excitation with zero gap, but this singlet 
is irrelevant. The transition at gcA is probably of first 
order, but is very close to second order. 

This work forms part of a research project supported 
by a grant from the Australian Research Council. We 
thank V. N. Kotov and J. M. J. van Leeuwen for very 
helpful stimulating discussions. We are also grateful to S. 
Sachdev for very important comments. The computation 
has been performed on Silicon Graphics Power Challenge 
and Convex machines. We thank the New South Wales 
Centre for Parallel Computing for facilities and assistance 
with the calculations. 



[2: 

[3: 
[4: 

[5] 
[6] 
[7] 

[s: 
[9: 
[10: 

[11 

[12 
[13 

[14 
[15 



S. Sachdev, in Low Dimensional Quantum Field Theo- 
ries for Condensed Matter Physicists, edited by L. Yu, S. 
Lundqvist, and G. Morandi, World Scientific, 1995. 
E. Dagotto and A. Moreo, Phys. Rev. Lett. 63, 2148 
(1989); F. Figueirido et al, Phys. Rev. B 41, 4619 (1989). 
H.J. Schulz, T.A.L. Ziman, and D. Poilblanc, J. Phys. 
(Paris) I 6, 675 (1996). 

N. Read and S. Sachdev, Phys. Rev. Lett. 62, 1694 
(1989); G. Murthy and S. Sachdev, Nucl. Phys. B 344, 
557 (1990). 

V. N. Kotov, J. Oitmaa, O. P. Sushkov, and Zheng Wei- 
hong, Phys. Rev. B 60, 14613 (1999). 
M.P. Gelfand, R.R.P. Singh, and D.A. Huse, Phys. Rev. 
B 40, 10801 (1989); M.P. Gelfand, ibtd. 42, 8206 (1990). 
R.R.P. Singh, Zheng Weihong, C. J. Hamer, and J. Oit- 
maa, Phys. Rev. B 60, 7278 (1999). 

E. Lieb, T. Schuhz, and D. Mattis, Ann. Phys. (N.Y.) 
16, 407 (1961). 

L. Capriotti and S. Sorella, Phys. Rev. Lett. 84, 
3173(2000). 

M. S. L. du Cro o de Jongh, J. M. J . van Leeuwen and 
W. van Saarloos, 



:ond-mat/0002116 



V. N. Kotov, M. E. Zhitomirsky, and O. P. Sushkov, 



cond-mat/0001282 



R. Guida and J. Zinn- Justin, J. Phys. A 31, 8103 (1998). 
A possibility gc2 < Qci has been excluded by the analysis 

i- 

V. N. Kotov and O. P. Sushkov, Phys. Rev. B 61, 11820 
(2000). 

J. Oitmaa and W.H. Zheng, Phys. Rev. B 54, 3022(1996). 



FIG. 1. Schematic picture of the simple columnar dimer- 
ized state. The ovals represent two spins coupled into a sin- 
glet. 



FIG. 2. The plot of 1/X-Pi where xp is the plaquette sus- 
ceptibility calculated in the simple columnar dimerized state 
using dimer series expansion. The value of 1/xp vanishes 
at gc3 = 0.50 ± 0.02 indicating a transition to the columnar 
dimerized state with plaquette type modulation. 



FIG. 3. The plot of 1/xd, where xd is the dimer suscep- 
tibility calculated in the simple Neel state using Ising series 
expansion. The value of 1/xd vanishes at gd — 0.34 ± 0.04 
indicating transition to the Neel state with spontaneous spin 
columnar dimerization. 



FIG. 4. Plot of the difference AC = 2|(Si • Sj) - 3(3^ S])] 
for 1st neighbors (full line), 2nd neighbors (short dashed line), 
and 3rd neighbors (long-dashed line) versus g = J2/J1. The 
full rotation symmetry of the ground state is restored when 
AC = 0. 
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FIG. 5. Schematic phase diagram and the excitation spec- 
tra of the Ji — J2 model. Solid lines show the triplet gap, the 
dashed lines show the gaps of the relevant singlets, and the 
dotted line shows the gap of the irrelevant singlet. 
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